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Abstract: 

We study a quantum analogue of the iterative perturbation theory by Kolmogorov used 
in the proof of the Kolmogorov- Arnold-Moser (KAM) theorem. The method is based on 
sequent canonical transformations with a "running" coupling constant A, A 2 , A 4 , etc. The 
proposed scheme, as its classical predecessor, is " superconvergent" in the sense that after 
the nth step, a theory is solved to the accuracy of order A 2 . It is shown that the 
Kolmogorov technique corresponds to an infinite resummation of the usual perturbative 
series. The corresponding expansion is convergent for the quantum anharmonic oscillator 
due to the fact that it turns out to be identical to the Pade series. The method is easily 
generalizable to many-dimensional cases. The Kolmogorov technique is further applied to 
a non-perturbative treatment of Yang-Mills quantum mechanics. A controllable expansion 
for the wave function near the origin is constructed. For large fields, we build an asymptotic 
adiabatic expansion in inverse powers of the field. This asymptotic solution contains 
arbitrary constants which are not fixed by the boundary conditions at infinity. To find 
them, we approximately match the two expansions in an intermediate region .We also 
discuss some analogies between this problem and the method of QCD sum rules. 



1 Introduction 



The remarkable progress made in a last few decades in the study of dynamical systems 
in classical mechanics has stimulated a renewed interest and revision of fundamentals of 
quantum mechanics. In particular, a number of efforts has been undertaken to elucidate 
possible manifestations of a (semi-) quantum chaos which, due to the correspondence prin- 
ciple, is expected to persist in this or that form in quantum mechanics, and to translate 
the concepts of non-linear classical mechanics to the quantum language ( see e.g. Q for a 
review ). One may hope, at the same time, that powerful calculation schemes of classical 
mechanics, being respectively transformed, could yield new methods supplementing more 
traditional tools of quantum physics such as the perturbation theory. 
One such an approach is the Kolmogorov super convergent method of classical mechanics 
|^] which has been the working instrument in the proof of the celebrated KAM theorem || 
. In a shortened form, the main Kolmogorov's idea can be introduced as follows. Imagine 
we have added some non-integrable perturbation eV\ (1, 6) to a integrable Hamiltonian 
Hq(I) which is the function of the action variables only. It has been known for a long time 
0] that successive canonical transformations (1,0) — > — > ... — > (I n ,9 n ) are able to 

lower the formal order of the perturbation by one order of the small parameter e at every 
step eV± — > e 2 V2 — ► ... — ► e n V n . However, this procedure generally results in a diverging 
asymptotic series applicable only at fairly short time intervals. The divergence originates 
in a possibility for an existence of rational ( and close to rational ) resonant relations 
among the frequencies uji of the unperturbed motion YU m i u i = with a set of the integer 
mi's . This is the so-called small denominators problem ( we refer the interested reader to 
the classical textbooks ||, Q and the excellent review paper |7| on the whole subject of 
the classical perturbation theory and the KAM theorem ). As any irrational number can 
be approximated to an arbitrary accuracy by rational numbers, one can readily imagine 
the scale of difficulties to be met with. It was Kolmogorov Q who has managed to show 
for the first time that successive canonical transformations may be chosen such that the 
order of the perturbation is increased by the square of the proceeding one for each step 
: eV\ — ► e 2 V2 — ► e^V^... — > e 2 " V n . This phenomenon of superconvergence occurs as 
a result of a clever redefinition of what is to be named "unperturbed" and "perturbat- 
ing" Hamiltonian at each iteration step. This idea actually dates back to the well known 
Newton method for finding the root of the equation f(x) = . Instead of expanding the 
function into the Taylor series around a suspected point xq , the Newton method suggests 
a simple first order iterative procedure of the form x n+ \ = x n — £fe4 • Thus, the value 

J \ x n) 

about which the expansion is made is moved closer to the true root at each iteration, that 
results in the superconvergence of the Newton method ( more on the analogy between the 
Newton method and the Kolmogorov approach can be found in § ). It turns out that 
owing to the superconvergence, the new perturbation theory is able to overcome the effect 
of near-to- resonance small denominators. Consequently, we end up with the convergent 
procedure, provided the initial conditions are adjusted at each step such that not to get 
just on a resonance surface || , f|] , || , @]. It has been shown || how to estimate from below 
the smallness of the resonant denominators that represent a set of frequencies of motion 
of a measure zero. The KAM theorem §,[§ states that most non-resonant invariant tori 
are only slightly deformed by small smooth Hamiltonian perturbations and " form the ma- 
jority in the sense that the measure of the complement of their union is small when the 
perturbation is small" ||. 
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In this paper we study a quantum analogue of the Kolmogorov technique. We find that, 
besides of the obvious aesthetic attractiveness (the superconvergence in the above sense) , 
the novel scheme possesses other interesting properties. In particular, it will be shown 
that the constructed super convergent procedure corresponds to an infinite resummation 
of the usual perturbative series and in some (non-trivial) cases yields a convergent quan- 
tum iterative scheme. Another observation is that the super convergent method can be of 
use in cases when the usual perturbation theory fails. To exemplify this property, we work 
out so-called Yang-Mills quantum mechanics ||]. Besides being a model for an infra-red 
behavior of the YM fields Q , this system is of particular interest due to the fact that 
its classical counterpart possesses strong chaotic properties |L0| . Note in this respect 
that the most popular nowadays approach to the quasiclassical quantization of classically 
chaotic systems, based on a summation of periodic orbit [|J , may be not the best way 
in view of its formidable complexity (the number of periodic orbits grows exponentially 
with the energy in a system with hard chaos). Sometimes it can be expected that the 
corresponding purely quantum system is much easier to solve. Then, given a quantum 
solution, the quasiclassical limit can be followed up. Finally, we mention that though our 
presentation is restricted by quantum mechanical examples, the formalism may hopefully 
be extended to the field theory. 

The presentation is organized as follows. In Sect. 2 we solve the standard textbook problem 
of building the quantum mechanical perturbation theory in a somewhat unusual manner 
which elucidates the quantum superconvergence phenomenon. The practical inconvenience 
of this approach turns out to be similar to that of the classical perturbation theory based 
on canonical transformations (see below). Then in Sect. 3 we present an equivalent but 
more convenient "non-coordinate" approach based on the so-called Lie transformation 
technique (see e.g Q ). In this form, the method admits ( at least, formally ) a generaliza- 
tion to many-dimensional cases. Its additional advantage is a formal equivalence between 
the classical and quantum perturbation theory. The Kolmogorov technique within the Lie 
transforms is further outlined. The proposed scheme is applied in Sect. 4 to the classical 
anharmonic oscillator and the reconstruction of the perturbative series is demonstrated. 
The corresponding problem for the quantum anharmonic oscillator is studied in Sect. 5 . 
The obtained series is shown to be equivalent to the well known Pade series, thus providing 
convergent expansions for both small and strong couplings. We check by direct calcula- 
tions that a re-expansion of the Kolmogorov series back into the standard power series 
reproduces the known results for the ground state energy. Finally, we apply in Sect. 6 the 
super convergent method to physically interesting YM quantum mechanics. Avoiding the 
small denominators to lowest orders is demonstrated. We show that the modified pertur- 
bation theory yields a controllable approximation to the behavior of the wave function 
near the origin. In a large fields region, we make use of another approximation and obtain 
an asymptotic adiabatic expansion in inverse powers of a field. This solution contains 
arbitrary constants which are to be determined from a matching with the small distance 
expansion. We make the approximate matching in two ways. The first one is essentially 
the least square matching over a "stability region" similar to the procedure used in the 
QCD sum rules method [11]. The second approach is a matching of the wave functions and 
their derivatives in a distribution sense (see below). Both methods yield similar results. 
The ground state energy is calculated with an accuracy of order 3.5 % while the accuracy 
of matching the wave functions is worse, of order 15 % . Sect. 7 contains a brief summary 
and several concluding remarks. 
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2 Super-convergent perturbation theory 



Let us suppose we know the solution to the unperturbed problem 

H ^ = £<M 0) (1) 

and we try to solve the perturbed problem 

(H + XW )% l = E n ^ n (2) 

We would like to solve this standard textbook's problem by a nonstandard method fol- 
lowing as close as possible the Kolmogorov's idea of super convergence (see Sect. 3 for the 
corresponding classical treatment). To this end, let us look for a solution in the form of a 
unitary transform 

* n = (e-^0 (1) )n = *P - A £ (n^\Si\m^)^ + . . . (3) 

where Si is an anti-hermitean operator Sf = —Si. We try to choose Si in such way that 
the resulting equation for will be in some sense simpler than the initial equation (2). 
We obtain 

e XSl (H + \W )e- XSl (i) W = E<f>W (4) 



or 



{H + X(W + [Si,H ]) + X 2 W 2 + A 3 ^ 3 + A 4 ^4 + 0(X 5 ))(p^ = E<f>W (5) 



where 



W 2 = [Si,W ] + ^[Si,[Si,H ]} 

W 3 = ^[Si,[S u W }} + ^[Si,[Si,[S u H ]}} (6) 
W A = ±[Si,[Si,[Si,Wo}}} + ^[Si,[S u [Si,[Si,H }}}} 

Let us require that the operator Wq + [Si , Hq] has only diagonal matrix elements in the 
basis of the eigenfunctions of the operator Hq : 

(n {0) \W + [Si,H }\m i0) ) = 0,n^m 

(n^\W + [Si,H ]\n^) = (n(°)|^ |n {0) ) (7) 

It follows from (7) that 

/ (0), o | (0)\ <n (0) |WoN (0) ) , m 
(n^>\Si\m^>) = ^ j^-,n^m (8) 

E n — E m 

If we omit the terms 0{\ 2 ) and higher, we thus obtain the usual formulas of the first order 
perturbation theory. Indeed, to this accuracy the exact solution to the equation (6) is just 
the unperturbed initial wave function : the equation 

Hi<p n = (H + \(Wo + [SuHoDfa = E n <p n (9) 
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has the solution 

4> n = 4°) , E n = 4°) + \(nW\W \nW) (10) 
and thus the solution to our problem to the given accuracy is 

^) = ^) + AV < ro( ° )|Wb|w(0)) ^) (11) 

Let us now try to make more precise calculations. The difference from the textbook's 
solution comes due to the fact that we retain the prescription (8). Then there is nothing 
to prevent one treating the Hamiltonian (9) as the new unperturbed Hamiltonian in the 
new perturbation problem 

(Hi + X 2 W 2 + X 3 W 3 + X 4 W A + 0(A 5 ))</> (1) = E<tF> (12) 

It is very important to notice that the order of the perturbation in (12) is A 2 ( and higher 
powers ). The term 0(A) hidden in the Hamiltonian H\ is no more the small parameter. 
This is in a sense " another A " as it is now a part of the exactly solvable Hamiltonian H\ . 
One may temporary call it e.g. p in order not to confuse it with the real small parameter 
A 2 . As we will see in a moment, this redefinition leads to an infinite resummation of the 
Rayleigh - Schrodinger perturbation series. 

We would like now to make again a canonical transformation such that the transformed 
operators W 2 , W 3 etc. become diagonal in the basis of the eigenstates of the Hamiltonian 
Hi. Obviously, this could be done with the choice <f>^ = exp(— X 2 S 2 )(f)^ > then the term 
W 2 is diagonalized. However, the term W 3 can be also diagonalized at the same step. To 
this end, we have to choose 

(1) = e -A2S 2 -A3Sa (2) {s + = _ ^ g+ = _ Ss) (13) 

Let us show that with this choice we arrive at independent equations for the operator 
S2, S3 which are as simple as eq. (7). ( an attempt to diagonalize also the term A 4 H4 
would break this property ) . We obtain the following equation for the wave function <f>^ : 

(Hi + X 2 (W 2 + [St, Hi}) + X 3 (W 3 + [S 3 , Hi}) 
+X 4 (W A + [S 2 ,W 2 } + ^[S 2 , [S 2 ,Hi}}) + 0(X 5 ))^ = (14) 

Proceeding analogously to (7) ( and taking into account (10) ), we obtain 

/( o)| Q | (0)\ (n^\W 2 \m(°)) (0) (0) (n(°)|H/ 3 M 0) ) , n 
<n< ^|5 2 |m^ >) = — - , < n t >|S 3 K j ) = f( i) f( i) > n + m ( 15 ) 

(The tacit assumption we have done here is that the new perturbation problem in not de- 
generate. If this is not the case, our procedure terminates, and the degenerate perturbation 
theory must be used. However, the degeneracy is not possible in the one-dimensional case, 
provided the operators Hi, W\ are self- adjoint.) 

These canonical transformations can be continued. For example, at the third step the 
terms up to A 7 will be diagonalized in respect to the eigenbasis of the Hamiltonian 

H 2 = Hi + X 2 W 2 + X 3 W 3 (16) 
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where W2,W^ are the operator in front of the powers A 2 , A 3 in (14). What we have ob- 
tained is just the quantum version of the super convergent procedure by Kolmogorov ( see 
Sect.l and Sect. 3 ). A few comments are now in order. 

(1) Transferring the averaged terms from the perturbation to the new unperturbed Hamil- 
tonian is the most non-trivial part of our scheme. This is the precise analogue of the New- 
ton method and the Kolmogorov approach in classical mechanics. The only difference is 
in the averaging procedure. We will see in the next sections that this quantum averaging 
can be done formally identical to the classical one. 

(2) The suggested method corresponds to the infinite resummation of the usual perturba- 
tive series. Let us illustrate this obvious observation on the example of the 0(A 4 ) results. 
We obtain from (13) , (14) : 

#=#+) .j;fcW! + ... (17) 

We would like to emphasize once again that while we expand the nominator to the given 
accuracy, retaining of all powers of A in the denominator is completely legitimate, as they 
are coming from the exactly solvable problem with the Hamiltonian H\. Let us also show 
that rational functions of the coupling constant appear in the formulas for the energy 
levels starting from the third step. The equation to be solved is 

(H 2 + X 4 W 4 + . . .)(t> {2) = E(f) {2) (18) 

Proceeding as previously, we substitute (j)^ = exp (-A 4 5 4 + . . .)^ (3) . Then the correc- 
tions to the energy of order 0(A 4 ) are expressed in terms of the matrix elements of the 
operator W4 ( see (14) ). Thus we immediately observe an appearance of A in denomi- 
nators owing to (15) . In effect, we get automatically a Pade-type series for the energy 
levels. It will be shown in Sect. 5 that the corresponding series is convergent for the quan- 
tum anharmonic oscillator. 

(3) It can happen that at some step our equation for the operators Si will not possess 
anti-hermitean solutions, but instead will give raise to some symmetric operators. Then a 
care must be taken to define self-adjoint extensions of operators at hand since otherwise 
the eigenvalues would not be preserved by a new canonical transformation. 

(4) The practical inconvenience of the method as it stands is related to the necessity to 
re-calculate at every step the matrix elements from one to another basis. This problem is 
quite similar to an analogous problem in the classical perturbation theory H where the 
generating function is a function of both the old and new variables while the new Hamilto- 
nian must be expressed as a function of the new variables only. An elegant solution to this 
problem consists in proceeding to a "non-coordinate" description of the canonical trans- 
formations which is called the Lie transformation method ( see Sect. 3 ). More important 
is the fact that in such a form, the method admits a straightforward generalization to an 
arbitrary dimensional case. 
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3 Superconvergence with Lie transforms 



The classical perturbation theory reduces the formal order of the perturbation by virtue 
of successive canonical transformations chosen such that to eliminate the phases (this 
procedure is to be done with preventing secular terms, see ||, ||). A generating function 
S(Ii,9,t) contains the old (9) and the new (Ji) variables, thus a relation between the 
new Hamiltonian H(I\,9\,t) and the old Hamiltonian H(I,9,t) appears also in a mixed 
representation : 

H(I 1 ,9 u t) = H(I,9) + dS{h d f> t) (19) 

where 

r 9S , „ dS 

'-ae""""^^ <20> 

Eq. (19) is just the relation between the values of two Hamiltonians at corresponding 
points in the phase space. In order to convert it into a relation between the functions, 
one has to solve the functional equations (20). This is not very convenient when one 
goes beyond the leading order. An elegant modified version has been proposed by Hori 
and Deprit (see f| , p2| and references therein). In this version the transformation from 
the old canonical variables x = (pi,qi) to the new canonical x = (pi,q~i) is given by a 
Lie generator (rather than generating function) w(x, A) = u>x + Xw2 + . • • such that the 
transformation x — > x(x, A) is just the shift by the "time" A along the trajectories of the 
system with the "Hamiltonian" w : 

I - M < 21 > 

(we use the same marking for the Poisson bracket in (21) as for commutators in Sect. 2 
because, as will be clear soon, the formal structure of the corresponding classical transfor- 
mation turns out to be identical to that of the quantum problem) . One further introduces 
the Lie operator L , simply related to the vector field X w associated with the Hamiltonian 
w : 

L =K , = » _ = (22) 

dq k opk opkOq k 

and the evolution operator T which transforms any function g at the new point x(x, A) 
into another function / at the original point x : 

/[x] = 5 [x(x,A)] ^ f = Tg; x = Tx (23) 

The operator T is exactly what is needed in order to convert (19) into a functional equation 
since (23) means, in particular, that formally 

H = T~ l H (24) 

(we are only interested here in time- independent problems. For non-autonomous systems, 
the form (24) is incorrect, see ||, |i~2]| ). The operator solution for T can be readily deduced 
from the above formulas ; from (21) and (23) we obtain 

dT f x 

=- = -TL ==> T = exp[- / L(\ )d\] (25) 
dX J 
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To see the formal similarity of the presented formulas with those of quantum canoni- 
cal transformation [ 13 1 , let us retain for simplicity only the leading term of the power 
expansion L = L\ + \Li + . . . in (24). Then (see (35) below) 

x = T _1 x = (1 + ALi + \>?L\ + ...) = e ALl x (26) 
The inverse transformation is 

x = Tx = (1 — \L\ + ^A 2 L? -...)* = e~ XLl x (27) 
The Hamiltonian equations of motion written via the Hamiltonian vector field Xh 

dX -X x- Y( 9H - dH - )x (28) 
dt H ~ k d 1kdpk dp k dq k 

now read (we suppose L to be independent of time) 

e- XLl ^ =X H e- XL '5c (29) 

or 

§ = X S * (30) 

CO yn 

X R = e XL ^X H e- XL ^ = J2-r[Ll,[L 1 ,...[L 1 ,X H ]}] (31) 

n=0 n! v ' 

n times 

(the brackets mean the commutation operation when stand with vector fields) 
Using the identity [X w ,Xh] = — X\ w m || , we obtain from (31) the transformed Hamil- 
tonian 

OO yn 

tf=£^KK.. >,#]]] (32) 

^ Q n\> ■ 

n times 

The formal identity with the quantum mechanical formulas becomes quite transparent. 
All formulas of this section will be true in the quantum case, provided one substitutes 
"Poisson brackets" — > "quantum commutators" and "functions on the phase space" — > 
"operators in the Hilbert space". 
To get the perturbative series, one further expands 

OO OO 

n=0 n=0 

OO OO 



where 



T = ^A"T n ; T-^Y^^Tn 1 (33) 

n=0 n=0 

OO OO 

h = x n H n ; s = \ n H n , 



n=0 n=0 

where 



L n = [w n , ] ; T = l ; T^ 1 = 1 (34) 
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Then a simple algebra leads to the recursion relation (see [6] ) 

n-l 

n 



n— 1 

= — / > L n ^ m T (35) 



and, as the consequence of (24), ( n > 



m=0 



n-l 



Ai^n = "-(#n - H n ) - (L n - m H m + mT n } m H m ) , (36) 

m=l 

where 

D = [ ,H ] and H = H (37) 
To the fourth order, the above equation yields 

D wi = Hx-Hx (38) 
D w 2 = 2{H 2 -H 2 )-L l {H 1 + H 1 ) (39) 

D w 3 = 3(H 3 -H 3 )-L 1 (H 2 + 2H 2 )-L 2 (H 1 + ±H 1 )-^L 2 1 H 1 (40) 

D w A = 4(ff 4 - ^4) - (^3 + + kil + k)lfi (41) 

3d 3 b 

- L 3 Hi - (L 2 + L\)H 2 — L\H 3 — 3L\H 3 — L 2 H 2 (42) 

The standard (von Zeipel's) perturbation theory is reproduced in this way as follows. 
Given Hi, one chooses H\ such that secularities in the r.h.s. of (38) are eliminated, and 
then finds w\. At second order, one substitutes the found w\ into (39) and looks for a w 2 
which eliminates secularities in the r.h.s. of (39), etc. 

As has been stressed in ||, the Lie transforms formalism is more convenient for elucidat- 
ing the superconvergence than the standard method of canonical transformations. The 
only difference from the above described line of reasoning is a different choice of the Lie 
generators wi. 

The principal rule is exactly parallel to the method used by us in Sect. 2 : successive new 

(i) 

Hamiltonians are obtained by successive Lie transforms while equations for different Wn 
and Wm { the superscript (i) will generally denote the number of the iteration ) must be 
independent of each other. At each step, chosen as many as possible independent u/^'k, 
all remaining 's are set equal to zero. We describe two steps of the Kolmogorov scheme 
following ||. 

At the first step, there is only one equation (38). To kill the secular term, we set 
Hi = < Hi > (< > here denotes averaging over the "fast" angle variables. In Sect. 
4 we will, however, follow another averaging procedure ). Then is determined from 

D w{ 1] = -{Hi} ; = , i > 1 (43) 

where { } denotes the oscillating part. All others H^ are given by (39 - 42) with the 
constraint (43). One finds 

= ^[w?,(hP+Hi)] (44) 
BP = -[^(tf^ + ^Vi])] (45) 
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One can show that the general recursion formula reads ( n > 2 ) 

s% = dh^'^ - i [w? ^ ])] (46) 

At the second step, we absorb the averaged (i.e. solvable) part A < Hi > into the zero 
order part of the new "old Hamiltonian" iJW : 

= +\<H 1 > 
h[ 1] = (47) 
= HP , i>l 

After the second transformation , we obtain a new Hamiltonian with = H { q ] . 

(2) 

As there is no first-order perturbation, we may choose w\ = 0. Thus, we obtain from 

(2) (2) 

(39 -42) two independent equations for w 2 and w 3 : 

D^w^ = 2(SP-hP) (48) 

= HS^-hP) (49) 

while 

= 4{BP - H^) - [ W ?\(H^ + H?>)] (50) 

where the new time derivative is calculated in respect to the zero order new "old" Hamil- 
tonian : 

D^ = [ ,H<p] (51) 

- r<2) — (2) 

The equation (48), (49) are solved simultaneously choosing H2 and to eliminate 

(2) (2) (2) 

secularities in w 2 and , as has been done in (43). All w\ , i > 3 are set equal to 

zero. Then the next terms in the transformed Hamiltonian are 

= HV + \[ W f\HW+ H V] (52) 



-1 it — _i_ 

= H^ + Wi^-rnH^+rnT^H^) (53) 

m=l 

At the third step, the new "old Hamiltonian" is built analogously to (47) : 
flf = + A 2 < > +A 3 < HP > 

H[ 2) = , i = 2 , 3 (54) 
HP = HP , i > 3 

This iterative procedure can be continued. It can be shown by induction || that at the 
nth step, the terms from i = 2 n to i = 2 n+1 — 1 are diagonalized. This completes the 
description of the Kolmogorov method within the Lie transforms formalism. 
It should be mentioned that, while being extremely important theoretically (for the proof 
of the KAM theorem), the Kolmogorov method did not get much practical applications. 
The only known to the author example has been given by Chirikov Q who has shown 
how to make the first Kolmogorov's step in the problem of the classical one-dimensional 
pendulum using the usual formalism of canonical transformations. But, as we have seen 
in Sect. 2 , really interesting phenomena within the Kolmogorov method start only from 
the second step. It is also worth noting that the notion of superconvergence has been used 
in a different content in [13] in treating the Birkhoff-Gustavson normal form. 



9 



4 Classical anharmonic oscillator 



For this one-dimensional classical problem, the use of the perturbation theory is nothing 
but a formal exercise. Indeed, the corresponding equation of motion can be solved in terms 
of elliptic functions that indicates that a perturbative series is convergent. Obviously, this 
convergence is related to the absence of resonances in any one-dimensional system. The 
perturbative series in this case is a formal expansion of the Hamiltonian in powers of the 
harmonic oscillator Hamiltonian. A final answer is usually expressed in the form of an 
expansion of the frequency of non-linear oscillations in powers of the energy. We will study 
this problem merely for the sake of comparison and working out reliable tools for treating 
the corresponding quantum problem. 
We start with the Hamiltonian 

H(p, x) = ^p 2 + ^uj 2 x 2 + \uj 2 x 4 (55) 

It is very convenient to proceed from the canonical variables (p, x) to the new canonical 
complex-conjugate variables 

1 



{ujx + ip) 

o 

z = , — (lux — ip) (56) 
s/2u 



Then 

H = # + A#i 

Hq = lozz (57) 

Hi = -(z + z) 4 = -(z 4 + Az 3 z + Qz 2 z 2 + Azz 3 + z 4 ) 
4 4 

In the complex variables the Poisson bracket takes the form 

r , i -fdfdg dfdg 

[f ' 9]= - t{ d-zdz-Wzd-z } (58) 

The advantage of proceeding to (z, z) and not to the action-angles variables is three-fold 
: (1) it makes quantum mechanical formulas written in the second quantized form very 
similar to classical ones; (2) simple explicit formulas for the Lie generators to all orders 
can be given ; and (3) the averaging procedure in the complex variables becomes much 
simpler in comparison with averaging in the action-angle variables. Instead of integration 
over the angles, the averaging is simply given by picking up terms with equal powers of 
z and z in a product. This procedure is obviously equivalent to the integration over the 
angle variables. 

Proceeding along the line of Sect. 3 , we obtain at the first step 

Hq 1 ^ = Hq = ujzz 

= <H l >=h 2 z 2 (59) 



Solving eq. (43) yields 



w\ 1] =i^-(\z 4 + 2z 3 z-2zz 3 -\z 4 ) (60) 
1 4uj 4 4 
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The transformed Hamiltonian is calculated according to (44), (45) (46) : 

= -L[-34z 3 z 3 + (z 6 -6z 5 z-33z 4 z 2 + (c. c. ))] (61) 
8u> 

= ^-^[2250z 4 z 4 + {33z 8 + 72z 7 z + 348z 6 z 2 + 1272z 5 z 3 + {c. c. ) )] (62) 

= 5 \-52Mz 5 z 5 + (37z w - Wz 9 z - A35z 8 z 2 - 17Q0z 7 z 3 

2oou> a 

- 4210zV + ( c. c. ) )] (63) 

where ( c. c. ) means the complex conjugation. The lengthy intermediate calculations have 
been made with the help of the Mathematica package of analytic calculations. Proceeding 
to the second step, the new "old" Hamiltonian becomes 

= cjzz + \-z 2 z 2 
H[ 1} = (64) 



We choose 



r(2) _ / o-(l) ^_ 17 „ :■ 



H-2 = <H?> >=-— (65) 

H® = < H i i ) >=—z 4 z 4 (66) 
3 3 16w 

The equations (48), (49) are of the form 

[uM 1 '] = u^fU + |^V] = -n{ffW} (67) 
To solve these functional equations, we use the following formula 

H§>] = iuj{m - n)z m z n (68) 



1 + ^zz 



(2) (2) 

which can by checked by the direct calculation. Then answers for w\ ,w\ read 



(2) • 1 1 A -6 3_ 5 33 _ 4 2 



-z b - -tz - -z 4 z 2 - ( c. c. )] (69) 



4w 2 l + 2Azz6 2 2 



(2) 



1 1 r 33 



48lu 3 1 + ^zz 1 8 



[—z 8 + I2z 7 z + 87z 6 z 2 + 636z 5 z 3 - ( c. c. )] (70) 



We give the final form of the Hamiltonian to the 0(X 7 ) accuracy ( Here I stands for the 
product zz ) : 

= w/ + ^ A / 2 -^-/ 3 + ^4 /4 
2 4 oj 16 or 

15 A 4 /5 3563 + 13496^/ + 11799^ I 2 ^ 



320 lo 3 (i + 3 A/)2 

1 A 5 .285299 + 949266^1 + 566055 4/ 2 
I I ■ u 

320 (l + 3^/) 2 

1 A 6 7 4239854 + 33800451^7 + 939122824^ 2 + 991621174^ 3 + 229513504i" 4 
512 u 5 ' ' " ' (i + 3 A/)4 
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Obviously, at the next steps we will get the successive zero order Hamiltonians in the form 
of rational functions of fjj . The corresponding u;W's can be calculated using the following 
generalization of (68) : 

Q(zz)] = -i(n - m )z m z n ^—± (72) 



L P(zz) ' ^ v n v ' P(zz) 

where P{y) , Q{y) are arbitrary functions, P(y) 7^ and the derivative Q'(y) exists. To 
solve the equation 

[w{z,z) , Q{zz)\ = p^~\ (73) 

appearing in higher orders, we have to choose in (72) 

P(zz) =Q'(zz)P 1 (zz) (74) 

The equations (72-74) in principle solve the problem of finding the corresponding Lie gen- 
erators v)® at any order. There are a few lessons to be learnt from the above calculation. 
In the second quantization representation, which is the relevant arena for building a corre- 
sponding quantum scheme, one can expect somewhat analogous to (71) expression for an 
approximate quantum Hamiltonian, the action / being substituted by the particle number 
operator N = a + a , i.e. a kind of the operator Pade expansion. The obtained Hamiltonian 
has to be checked for being self-adjoint. Furthermore, to proceed beyond the first step, 
quantum substitutes of eq. (68), (72) are needed. 



5 Quantum anharmonic oscillator 



The quantum anharmonic oscillator is the standard test problem for any new quantum 
field theory - oriented method. A main feature, that the anharmonic oscillator shares with 
field - theoretical problems, is the zero radius of convergence of the Rayleigh-Schro dinger 



perturbation series [14] which prevents an extension of perturbative results into the strong 



coupling region. There has been a number of attempts to build strong coupling expansions 



meant to be generalizable to the field theory, refs. IS] may probably be mentioned among 
the most interesting ones. 

Intuitively, the divergence of the perturbative series can be understood as a result of the 
inability to reproduce a behavior in the vicinity of the essential singularity at the origin 
( the Dyson singularity ) in terms of an expansion in powers of the coupling constant. 
However, it is well known that the convergence can be gained after the perturbative series 



is resumed, the most popular resummation methods are due to Pade and Borel, see e.g. [16] 
for a review. In this section we would like to argue that the quantum superconvergence 
leads to the uniform convergence of the corresponding series for both small and strong 
couplings due to the fact that our results have the form of the (operator and functional) 
Pade series. 

For the quantum version of the Hamiltonian (55) (with u) = 1) we introduce the creation 
and annihilation operators : 

(x + ip) , a = — (x — ip) (75) 



2h V2H 
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and the particle number operator N = aa satisfying the commutation relations 

[a , a n ] = na" 1 - 1 , [a n , a] = na 11 ' 1 (76) 
The unperturbed and perturbating Hamiltonians are 

H = h(aa + ^) (77) 

h 2 

Hi = — [3(l + 4aa + 2a 2 a 2 ) + (a 4 + 4a 3 a + 6a 2 + (/i.c.))] 



we obtain h[^ =< Hi > and then 



where (/i.e.) means the hermite conjugation. Proceeding analogously to the classical case, 

a 

w <f> = ^(Ia 4 + 2a 3 a + 3a 2 - {h.c.)) (78) 
The transformed Hamiltonian , according to (44-46), is 

H^ ] = ^ [-42 - 288aa - 306a 2 a 2 - 68a 3 a 3 + . . .] 

fiW = — [333 + 3582aa + 6291a 2 a 2 + 3000a 3 a 3 + 375a 4 a 4 + . . .] (79) 

h{ 1] = -^-[-26649 - 390366aa - 977562a 2 a 2 - 741636a 3 a 3 - 196650a 4 a 4 - 15732a 5 a 5 + . . .] 
256 

where we have explicitly shown only the diagonal parts of the operators, the complete 
forms are rather lengthy. Proceeding to the second step, we add the diagonal part of 
the perturbation Hi to the unperturbed Hamiltonian Hq to get the new "unperturbed 
Hamiltonian", just as we have done in the classical treatment : 

= h(N + ±) + ^f(l + 2N + 2N 2 ) 
h[ 1] = (80) 
H\ l) = BP , *>1 



At the second step, the quantum analogue of eq. (67) can be written as 



[N + kN 2 , W ?] = 2 -^l{H^ } (81) 



where we have denoted 

k = 2 „ 82 

i + ^xn y ' 

The quantum version of eq. (68) reads 

[N + kN 2 , 1 -———a n a m ] = (n - m)a n a m (83) 

L 1 - k(n - m) + 2kN x ' 

(2) 

The above formula yields the following answer for the Lie generator w\ 

(2) (l-k)h 2 A 1 _ 6 15 1 _ 4 171 1 _ 2 

«>2 = 5 \-~ ~- - ~ a - —- ~ - ~ a — 



l 3 1-6k + 2kN 2 1 - 4k + 2kN 2 l-2k + 2kN 

1:52 ra\~- *— a*a- -, j^ aV - (h,.)] (84) 



2k + 2kN 1 — 4Jfe + 2kN 1 - 2k + 2kN 
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(2) 

The result for the Lie generator u> 3 has a similar structure. To calculate higher order 
terms of the transformed Hamiltonian, one needs a few more commutators which can be 
obtained as the simple consequences of the Haussdorff identity and Schwinger exponential 
representation of operators ( a and f3 are arbitrary c-numbers ) : 

[ ^ a ' m ' ^TPN ] = ( a-/3(n-m)+/3iV " ^TW^^ (85) 



TO— 1 | 

E 777 tttt— ^aV (86) 

~* , Z!(m-Z)!(n-m + Z)! 

!=mai{0,m-n} 



(2) 

Using these commutators, one can easily check that the operator w\ is indeed anti- 
hermitean. As the resulting expression for the second- iteration Hamiltonian is rather 
lengthy, we will retain below only the terms possessing non-zero vacuum expectation 
values : 

ffW = - + -Aft 2 - ^A 2 ft 3 + ^A 3 ft 4 - ^A 4 ft 5 
2 4 8 16 256 

A 4 ft 5 (l - k) r 29241 1 ^675 1 i 30 1 



16 L 16 l + 2k + 2kN 2 1 + 4k + 2kN l + 6k + 2kN> 

444609 A 5 ft 6 (l - k) 208791 1 

+ A 5 ft 6 H (87) 

512 4 L 32 l + 2A; + 2A;Af v ; 

14985 1 315 

+ 8 1 + 4k + 2kN ~ 1 + 6k + 2/fcA^ + " ' 

One can easily see that the perturbative expansion for the ground state energy [fbf l 

ft 3 xt2 21 - . 333 . 3 . 4 30885 x4t5 916731 5t6 

E = - + -A/i 2 A 2 ^ H A 3 ?i 4 A 4 ft 5 H A 5 ft b + . . . (88) 

2 4 8 16 128 256 v ; 

is reproduced from our formula by re-expanding the rational fractions. For example, 
about 25 % of the perturbative result at the A 5 order is due to the expansion of our 0(A 4 ) 
fraction. It seems plausible to suggest that the Kolmogorov series reproduces the usual 
perturbative result at all orders, though we do not get a general proof. 
Heuristically, one can expect the convergence of the Kolmogorov series (87) for the ground 
state energy due to the fact that, the higher is the order of the perturbative coefficient, 
the larger number of fractions will contribute to its value. This statement can be given a 
formal proof. The resulting from our procedure fractions can be reorganized into rational 
functions of the coupling constant. A convenient choice is e.g. to collect into a single 
rational function all fractions resulting from one Kolmogorov's step. Let us suppose that 
at some step we have obtained the answer in the form 

JVii(A) = f{Ni,Ni+ji](^X) = R[Ni+M /q[Nx\ (89) 

where B.W(q[M]) 

is a polynomial of degree N (M). (note that in our formulas j\ > 0) and 
at the next step the answer is of the form 

P N2 , j2 (X) = f [N2 ' N2+h] (X) , (N 2 >N 1 ) (90) 

Let us prove that j\ = ji. Suppose this were not true. As PN x j\ > Pn 2 ,j2 can be- thought 
as the representatives of two series of the Pade approximants, there exist the limits pX 



/ ix (A)= hm , / j2 (A) = lim fl N > N +^ (91) 
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and the convergence is uniform in the cut A plane. However, fj ± = fj 2 since a mea- 
sure giving the perturbative coefficients in the form of its moments is unique for the 
anharmonic oscillator problem 1C]. Then it means that the series fl N ' N+ j^ — fl N > N +h] 



converges uniformly to 0. As j\ , ji > , this is only possible when j\ = ji (at very large 
A fl N < N +fi ~ \i ). The fact that ji = j2 means that the Pade series obtained from Kol- 
mogorov procedure is diagonal and therefore converges uniformly to the true eigenvalue 
p7| . (Note that we have explicitly checked that all poles resulting from the Kolmogorov 
series are located at real negative A for any energy level and the resulting Hamiltonian is 
self-adjoint.) 

We would like to end up this section with a few comments. At any order, the eigenfunc- 
tions of a resulting effective Hamiltonian are the harmonic oscillator wave functions in 
proper coordinates. The proposed approach is essentially the method of a free field rep- 
resentation. Given the needed accuracy, the interaction is transferred to such canonical 
transformation that yields the free (i.e. harmonic oscillator) Hamiltonian. The perturbed 
wave function in the initial coordinates is given by the inverse transformation of the basis. 
The theory of linear canonical transformations in the Bargmann-Fock space is well studied 
[18]. As we deal with non- linear canonical transformations, more sophisticated methods 
are needed. We feel that a method suggested recently IS] can be of use here. Note also 
that the corresponding treatment in a many-dimensional case provides a non-linear exten- 
sion of the famous Bogolyubov transformations. 

The results obtained so far can be summarized as follows. In all cases where the usual per- 
turbation theory is applicable, the superconvergent methods corresponds to a Pade-type 
resummation of the perturbative series. In the above studied case, the convergence of the 
Kolmogorov expansion is the consequence of general theorems on the Pade approximation 
[16]. In other problems, convergence properties have to be tested separately. Instead of 
discussing this topic, in the next section we will show how a similar re-organization of the 
perturbation theory allows one to treat approximately a more complicated system with 
two degrees of freedom, where the usual perturbation theory cannot be applied. 



6 Yang-Mills quantum mechanics 



Yang-Mills (YM) classical mechanics has been proposed jl(J as a toy model for a study of 
non-perturbative nonlinear YM waves in the Minkowsky space. There one starts with fields 
in the SU(2) YM theory which depend only on the time (this is the reasonable assumption 
in the infra-red). Then, after the gauge fixing Aq = , the classical YM theory turns out 
to be essentially equivalent to the dynamical system with the Hamiltonian ficfl 




where the dynamical variables f a parametrize the gauge fields Af ( Of is a time-independent 
orthogonal matrix ): 

At = 0?f a (t)/g (no sum over a) and OfO- = 5 ab (93) 
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A simplified system obtained from (92) by fixing / 3 = is still highly non-trivial and 
commonly called YM classical mechanics: 

Hym = \{pI+P 2 2) + \x\xI (94) 

(we have introduced here the coupling constant A instead of 1/2 in (92) ). 

In spite of its seeming simplicity, the Hamiltonian (94) leads to a very complicated chaotic 

motion |1Q1 . It has been believed for a long time that the system (94) is ergodic [10]. More 



recently, small stability islands have been found p0| . There also has been a considerable 
interest in a quantum system with the Hamiltonian (94) [||, M, [21|, p2|] , ]23| (note that 



(94) is not the homogeneous-space model of the quantum YM theory, but still is closely 
related to the latter [||). The fact that a spectrum is discrete has been proved in pi] . 
Different methods, based on the quasiclassical and quantum adiabatic approximations |22f| 
and summation of adiabatically stable periodic orbits p3[] , have been proposed for finding 
the spectrum and wave functions. 

It is obvious beforehand that the usual perturbation theory is hopelessly inapplicable 
for treating the system (94). No matter how A is small, the perturbation results in a 
reconstruction of the spectrum from continuous to the discrete one, the huge effect. On 
contrary, a modified perturbation scheme a la Kolmogorov can be proposed which allows 
one to build a controllable approximation to the WF near the origin, where the adiabaticity 
breaks down (see below). 

Let us consider the "perturbation" in the second quantization representation (75) (we set 
h = 1 and use zt , Zi to denote the quantum annihilation and creation operators) : 

Xxjx 2 2 = \^[l + 2N 1 + 2N 2 + 4N l N 2 ] + \^{zjz 2 1 + zfz 2 i 

+ {l + 2N 2 )z 2 l + {l + 2N 1 )zl + {h.c.)} (95) 

In obtaining (95), we have used the non-commutativity of the quantum operators, that 
has allowed one to extract from the "perturbation" the part which is actually not a per- 
turbation as it is diagonal in the oscillator basis. The expression denoted by [...] in (95) 
must be treated as a part of an unperturbed Hamiltonian, if we intend to end up with a 
discrete spectrum. The fact that the "unperturbed Hamiltonian" is of the same order A as 
the remaining "perturbation" is not an obstacle since after the first step the perturbation 
is only of order A 2 (plus a small part of order A ). An additional problem arises due to 
the fact that both the unperturbed Hamiltonian and perturbation are degenerate so that 
the degenerate perturbation theory also cannot be used. We propose the following simple 
trick to handle this problem : let us just make both the unperturbed Hamiltonian and 
perturbation non-degenerate. We choose to transfer into the unperturbed Hamiltonian 
the terms 1 + 2N\ + AN\N 2 while the term 2N 2 is retained as the part of the perturbation. 
The concrete choice N\ «-> N 2 does not matter as the difference may be thought as the one 
to be repaired by higher order corrections. In other words, we sacrifice the exact symmetry 
of the problem in favor of only approximate one. (It is known that such procedure can 
yield wrong results for the large distance behavior, but, as will be shown below, in this 
region a different approximation becomes operative). As will be seen in a moment, this 
suggestion indeed leads to avoiding small denominators at the lowest order. Moreover, in 
such form the Hamiltonian is suitable for working out the Kolmogorov scheme, as the new 
perturbation contains the part which has a non-zero average (the 2N 2 term). 
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Adding to (95) the kinetic term operators, we arrive at the Hamiltonian 
H = Ho + eHi 

H = cJiiVi + uj 2 N 2 + /j,N 1 N 2 + a (96) 

#1=2^ ( 2iV2 + Z l Z 2 fnm + + Jnm^l ^2 + z l 9nmZ 2 + Z\ 9 nm Z 2 ) 



n,m=0 



where 

1 A 1 

uJi = - + -,uJ 2 = -,fi = X 

a = \ + \ ' £ = I ' 7 = 1 4 (97) 

and the only non-zero f'k , g'k are 

/ 02 = 7 + 2iVi , f 20 = J + 2N 2 , f 22 = l , 522 = 1 (98) 

(Note that the Hamiltonian of so-called Yang-Mills-Higgs quantum mechanics ||, which 
differs from (94) by adding the mass term ^uj 2 (xi+x 2 ) , can be treated along with the same 
line with a proper redefinition of parameters. The interest in the corresponding classical 
system has been recently renewed in connection with a study of non-abelian collective 
excitations of the quark-gluon plasma p4|.) 

To find the first step Lie generator , it is convenient to look for a solution in a form 



analogous to (96) [g5| : 

>? = E (^"ffi. - i^lT^T + *TcW zf - z?(c£lTz?) (99) 



n,m=0 



We further substitute (96), (99) into eq. (43) and equate the like powers of the quantum 
operators. The comparison of the coefficients provides 

a nm = fnm[nuJi + mto 2 + fi(nN 2 + mNi + nm)]^ 1 (100) 
c nm = 9nm[nui - muj 2 + fj,(nN 2 - miVi)] -1 (101) 

Note that the denominator in (101) never goes to zero. For the only nonvanishing g 22 = 1 
we obtain 

eg = [A + 2A(iV 2 - iVi)]" 1 (102) 

which does not turn to infinity at any eigenvalues of the operators N\ , N 2 . Using the 
machinery of the proceeding section, it is not difficult to calculate the next term in the 
approximate quantum Hamiltonian (again, we retain only the terms yielding non-zero 
vacuum expectation values) : 



= cj 1 (N 1 +N 2 )+fxN 1 N 2 + a + 



£ 



2 



(TV 2 - JVxXJVf - N 2 ) 



2 L wi + uj 2 - 4/i + /x(JVi + N 2 ) 
4 + 6(JVi + A^ 2 ) + 9A^ 2 + 2(N± + Af) + 3^1^ (iVi + N 2 ) + iVfjy| 

ui + w 2 + 2/x + /i(iVi + jV 2 ) 
(iV 1 2 -iV 1 )(iV| + 3JV 2 + 2) _ (N 2 2 -N 2 )(Nf + 3N 1 +2) 3 
wi - o; 2 + 2// + /i(iV 2 - JVi) ux - lo 2 - 2fi + fi{N 2 - iVi) J 1 M J 
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One can check by the term-by-term inspection that zeros do not appear in the denomina- 
tors in eq. (103). It can even be further speculated, similarly to the line of reasoning of 
Ref. [25], that small denominators do not show up at any order of the iterative procedure. 



For the ground state energy we obtain from (103) 

1,. A A 

1 — L _ 



£ ° = 2 (1 + 2-TTTU ) + 0(£ ) ' (104) 



or, setting A = \ 



£0=4(1 + ^-^)^0.611 (105) 



2 - 



This number is to be compared with the result of the numerical calculation |22| which 

reads Eq ~ 0.590. One sees that our accuracy is rather good, of order 3.5 %. 

Let us now turn to the study of the approximate wave functions (WF) (here we restrict 

ourselves by the first order formulas). The unperturbed Hamiltonian Hq (96) has the 

eigenfunctions 

= |nin 2 ) , s = (ni,n 2 ) , (106) 



or, in the coordinate space 



2 , 2 



where H n (z) stands for the Hermite polynomial. After a simple algebra we obtain the 
first order WF 



^(0) , 1 r Vni(ni-l)n 2 (n 2 -l) Ao) 

(rn.no) ^( n-i.no) o I 



(ni,n 2 ) ^(m,n 2 )- r 2 i UJl+u;2+ ^ ni+n2 _2y(ni-2,n 2 -2) 

+ 



y/ni(ni - l)(ra 2 + l)(ra 2 + 2) ( ) 



+ 



+ 



UJl - UJ 2 + H(n 2 ~ Til + 2) Y (ni-2,n 2 +2) 

^/(ni + l)(ni + 2)n 2 (n 2 - 1) (0 ) 

uj x - lo 2 + ii(n 2 - ni - 2) %i+2,"2-2) 
V(ni + + 2)(n 2 + l)(n 2 + 2) (0) 

Ul + 002 + H(m + n 2 + 2) <P (m+2,n 2 +2) V ) 

(7 + 2n 2 )^n 1 (ni - 1) ( ) _ (7 + 2n 2 )^(n 1 + l)(n x + 2) (0) 

wi+/in 2 * , (m-2,n 2 ) Wi+/in 2 %i+2,n 2 ) 



(7 + 2ni)y/n 2 (n 2 - 1) ,(0) (7 + 2m)y / (n 2 + l)(n 2 + 2) (p) , 

CJ 2 + ^ni ^(m.na-2) W 2 + jUTli < ? ) (n 1 ,n 2 +2)J 



The following from (108) asymmetry of the ground state WF (GSWF) is evidently an 
artifact of the asymmetry chosen in constructing the perturbation theory. As the GSWF 
must possess a highest possible symmetry, the final expression must be symmetrized in 
respect to the replacement (N\ «-> N 2 ) : 



+ (ii±M^M^) _ 27 ^ )1 + (ll „ I2)1 m 
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The obtained formula demonstrates that our expansion is controllable in a vicinity of the 
origin. With moving away from the origin, the correction constitutes about 50 % of the 
unperturbed value at x\ ~ 2.0 , if one moves along the lone X2 = , or x\ ~ X2 ~ 3.0 for 
the radial escape. ( Note that the validity of the ansatz for the wave function for small 
x, y in the form exp ^polynomial can be checked by the direct substitution of such an 
ansatz into eq. (110) ). 

To get approximate WF's in a large fields region, we need another approximation. To 
this end, we make use of an intrinsically true idea that, when one coordinate is large, 
the variables in the Schrodinger equation get adiabatically separated |22[ . Note, however, 
that the method used in Ref. [221 does not allow one to estimate corrections to the 



adiabatic separability, which is the question of great importance if we are going to match 
approximately two different expansions. 

Let us write down the Schrodinger equation (in what follows, we set x\ = x , X2 = y) 



1 d 2 d 2 

2 + dy^^^' y ) + Aa; V*(^ V) = E^(x, y) 



(110) 



and make in the region x > the point canonical transformation to the new adiabatic 
variables ( we are indebted to L. Franktfurt for this suggestion) : 



x 

y 



U = X 

t = ^y 



(111) 



The Schrodinger equation then transforms into 

l r d 2 ,d 2 2A , rix . . A t d It d 2 , , . 



where r) = n~4 ^(n, tu~2 ) ; the multiplier is due to the Jacobian of the transforma- 
tion (111). Note that in eq. (112) the variables get adiabatically separated in the region 
of very large u provided 



1 r 5$ 

4 u 2 dr 

t d 2 $ 
u du dr 





d 2 <& 




U dr 2 




d 2 <S> 


< 


U dr 2 



(113) 



In the asymptotic region we neglect in the zero approximation the r.h.s. and look for a 
solution in the form 



$o = £a r >M( u ) (114) 

ra=0 

where £ n (,T) satisfies 

[\^2 ~ Ar 2 )en(r) = -enUr) , (H5) 
i.e. the r - motion is the usual harmonic oscillator. The resulting equation for rj n {u) is 



-(2n + l)uun + 2Er] = 

du z 



(116) 
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where u> = v2A . The only physically accepted solution (damped as u — » oo) is given in 
terms of the Airy function |2q| : 



Tjn( u ) = const Ai 



2 F 1 

((2n + l)u;) 1/3 ('u- 



(2n + 1)uj' 

However, we are only allowed to retain the asymptotics of this formula, thus 



(117) 



$ (u, t) ~ A n e-%T 2 H n (T^j)^Te-i a ^-^) 3/2 , u -, oo (118) 

n=0 Ui 

where A n are the normalization coefficients and 



/ 2E 

a B = ^/(2„ + l) W , Pn = WT ^ (H9) 

Let us now come back to eq. (112) and try to calculate corrections. There is no room 
in (112) for the standard Green functions method because of the non-separability of the 
variables at arbitrary u,t. To proceed, we substitute into (112) the ansatz 

*(«, r) = $ (it, r)n(u, t) (120) 

and try to find the form of D,(u, r) in a vicinity of it — * oo. The boundary conditions on 
U(u, r) are 

f2(u, t)— > 1 , it — ► oo , r fixed (121) 
|fi(«, r)| < e^ r2 , t — ► oo , it /zxed (122) 



The equation for Q(u,t) reads (in what follows we denote $ = <&o , &t = etc.) 

(-*«r - 7^0^ + (2$ u + -$ r )ft u + -$fi UT 

u iu z u u 

+(2$ r u + - -^$)0 T + $n uii + $un TT = (123) 
u 4u 2 

According to the boundary conditions, we look for a solution in the form 

fl(u, r) = 1 + u-V^M + u-V2(r) + ^ 3/2 /3(r) + . . . (124) 

We further substitute this expansion into (122) and equate the like powers of u. The 
resulting equations for /i , /2 , fa read 

K(r)f[(r)+Ur)f^r)=0 (125) 
K(T)fL(T)+Mr)fZ(T) = (126) 
-a n <(r) + 2&(r)/£(r) + £ n (r)/£(r) = (127) 

Let us first concentrate on the equations (125), (126). It is easy to see that the non-trivial 
r - -dependent solutions are excluded due to the boundary condition (122). However, arbi- 
trary constants are not fixed by the boundary condition at the infinity. ( The acceptable so- 
lution for eq. (127) is a bit more complicated: fz{r) = C^ + a n J r dt^~ 2 (t) J * dzz£' n (z)£ n (z) 
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) . This means that our method is only giving a large-u parametrization of the wave func- 
tion. Recalling the Jacobian and returning to the original variables, we finally arrive at 
the GSWF of the form (we retain here only the n = term from (118) ) 

9{x, y) = Ari^-K-A)'" (1 + % + ^ + 2^ + • • •) (128) 

where A, C\, C%, C3 are some constants and II(r, C3) = C3 + \ J T dte Q o* 2 r(§, alt 2 ), T(j, z) 
stands for the incomplete Gamma function (we will omit this term in what follows) . 
Note that the obtained WF is in no sense an analytic continuation of (109). This is in 
the drastic contrast with one-dimensional quantum mechanics where the mathematically 



rigorous matching of the WKB wave functions can be given [26]. We would like to mention 



that the appearance of arbitrary constants in asymptotic large distance solutions is well 



known in hydrodynamics [27]. It is also known that small- and large-distance expansions 
can be analytically matched only if the small distance expansion has the infinite radius 
of convergence []2g|[ , which is not the case in our problem. Therefore, one has to rely on 
an approximate (or numerical) matching. In this respect, the situation is reminiscent of 
what one encounters in the QCD sum rules method (provided we discuss the matching 
problem for wave functions and not for Green functions). There two different asymptotic 
expansions, each of them being defined in its own range, are numerically matched in 



some intermediate region [11]. It may be useful to continue a bit this formal analogy. Our 



asymptotic adiabatic expansion is a method of parametrization of the large distance (fields) 



behavior and thus is analogous to the expansion over hadron states [11], the parameters 
A , G\ , C2 being the "hadron masses". The non-perturbative Kolmogorov technique 
substitutes the operator product expansion at small distances. Probably a bit deeper 
analogy with quantum chromo dynamics lies in the fact that the non-integrability of the 
system (94) blocks an explicit construction of " quasiparticle" weekly interacting operators 
in the intermediate region. It is not yet, however, clear whether the non-integrability or 
just the mathematical complexity of real QCD prevents one building hadron operators in 
terms of the fundamental quark and gluon fields. 

Inspired by this analogy, we have tried to find the coefficients A, C±, C2 by the numerical 
matching of the form (128) with the WF (109) at moderate x ~ 1 — 2 in the vicinity of the 
line y = . Though there exists a variety of matchings in three dimensions, the simplest 
method one can suggest is a minimization of the difference between the two wave function 
over the usual normalization measure with the integration done over the "stability slice 
" , i.e. a continuous limit of the usual least square method. As such procedure can yields 
jumps of the derivatives, it must actually be supplemented by some smoothing operation 
on the obtained functions. Within the least square method, this problem is just ignored, 
which is partially justified by the fact that we make the matching in the exponentially 
suppressed region where derivatives are of order of the functions themselves. 
We have found that there exists the "stability plateau" ranged from x ~ 1.5 to x ~ 2.5. 
With the matching in the intervals of x = [1.5, 2] or x = [2, 2.5] , the coefficients A, C%, C2 
change not more than by 10 % . The found coefficients are 

A = 2.50 ± 0.15 , Ci = -1.55 ± 0.10 , C 2 = 0.90 ± 0.10 (129) 

with the measure parameter for the interval x = [1.5, 2] 

f™dx f dy(AV) 2 M , . 

m = Jl - 20 ; ;) 2 - n% 130 

f^dxfdy® 2 
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where stands for the difference of the WF's (109) and (128). The measure parameter 
over the interval x = [2,2.5] is a bit worse, of order 13 % . The main defect of the 
above matching is the fact that it deals solely with the adjustment of the wave function 
without imposing any condition on the derivatives. As the wave functions are known only 
approximately, one cannot impose the condition of the continuity on the derivatives of the 
wave functions. However, in view of the approximate character of our solutions, in the 
intermediate region it seems reasonable to proceed to the description of the wave functions 
in terms of the distributions rather than the usual functions. This means that we introduce 
some smoothing which characterizes the actual accuracy of our calculations. Instead of 
the function at a point xq, we get the "function" at some small region Q containing the 
point xq : 

/(fi) = / f{x)<p{x)dx , (131) 

where 4>(x) is a differentiate function which is very strongly peaked at xo and vanishes 
outside the region fi. Then a precise meaning to the notion of matching two distributions 
fi and fi can be given. The distributions fi, coincide in a open set fi when 

/ M= I h$ (132) 

for every infinitely differ entiable function <j) whose support lies in fi pG| . The analogous 
conditions hold for the derivatives of /i,/2 : {fx,<j>') = if 2,4'') ■ F° r the purpose of 
the numerical matching, the test functions 4> should be strongly peaked at the stability 
plateau x = [1.5,2.5] whose existence is indicated by the least square matching. In our 
opinion, such matching in the distribution sense is adequate to the approximate character 
of both solutions. Note that a more ambitious way of solving the Schrodinger equation in 
the intermediate region would be a search for a test function eft satisfying the differential 
equation, i.e. a solution in the form of the integral transformation. 

For the approximate matching in the distribution sense, we have chosen the test function 
in the form of the smoothed 5 function 

cj)(x, y) = exp[— n 2 (x — xq) 2 — m 2 y 2 ] (133) 

7T 

and varied the parameters n, m, xq in the vicinity of the values n = m = 1 and xq = 2. 
Then the parameters A, C\, C2 have been determined from the matching of the integrated 
wave functions and the their partial derivatives over x ( the matching of the partial deriva- 
tives over y is the triviality when one chooses the test function (133) and neglects the term 
n(r, C3) in (128) ). As in this case there are only two equations for three unknown pa- 
rameters, we have taken one of them as the known one from (129) and found the other 
ones. We have found the same values of the parameters with the accuracy quoted in (129). 
We believe that the consistency between the two methods of the matching confirms the 
reliability of the results (129). Still, as we have seen, these two methods are not quite in- 
dependent. Another method of a matching could be a search for global conservation laws, 
that would allow one to find the parameters A, C±,C2 or some relations among them. 
Unfortunately, we have not succeeded in this attempt. 
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7 Concluding Remarks 



In this paper we have studied a quantum extension of the Kolmogorov iterative technique. 
We have found explicit solutions to functional equations determining the high orders Lie 
generators and Hamiltonians in both the classical and quantum cases and tested the whole 
approach on a few problems. For YM quantum mechanics, we saw that the Kolmogorov 
technique for the small distance regime must be supplemented by the adiabatic large dis- 
tance approximation, thus the matching problem arises. As the whole affair is still in its 
infancy, one can hardly say beforehand whether it is able to suggest a reliable calculation 
scheme for more realistic problems of quantum mechanics and quantum field theory. Nev- 
ertheless, as we have seen, the KAM technique possesses a few very attractive features. 
The success in the construction of the convergent iterative scheme for the quantum anhar- 
monic oscillator may allow one to hope that this technique could also be helpful in more 
complicated problems, provided the method is easily generalizable to arbitrary dimensional 
cases. The hope to overcome the Dyson instability is due to the non-linear character of 
the KAM technique. Moreover, because of its recurrent definition, this method seems 
to be ideally suited for an analytical computerization. The need in the computer's help 
is seen already at first Kolmogorov's steps. We want to emphasize, however, that these 
difficulties are only technical and may be easily handled. On the other hand, the KAM 
technique is essentially non-perturbative and suggests the iterative non-linear extension 
of the Bogolyubov transformations method, when applied to infinite dimensional systems. 
Thus, it can be of some interest for clarifying a connection between perturbative and 
non-perturbative effects in quantum field theory. 



The proposed technique could be further developed in a few possible directions. One 
of them is working out a theory of non-linear canonical transformations [ 19 1 , which would 
allow one to calculate approximate wave functions directly from the "non-coordinate" Lie 
transforms formalism. Another problem is an improvement of the matching procedure 
and the use of a more accurate smoothing technique analogous to that used in classical 
mechanics. Also a search for global conservation laws seems to be rather important for the 
matching problem. The matching procedure done in this paper is very crude and, in par- 
ticular, leads to an admixture of excited states to the ground state wave function. Surely 
enough, the matching must be improved to make this method work. Such an improvement 
would be especially important for the matching large- and small- distance solutions for 
highly excited states, where chaotic properties are expected to manifest themselves in full. 
It would also be rather interesting to reformulate the present approach in the path integral 
formalism and apply it for a calculation of Green functions. Hopefully, this could be done 
in the coherent states representation of the path integral. In such a form, this method 
could be directly confronted with the Operator Product Expansion, that might lead to a 
deeper understanding the phenomenological success of the QCD sum rules method [11]. 
As for the problem of calculation of the spectrum, the generalization of the present ap- 
proach is quite straightforward. In particular, we hope that this method may be extended 
to calculate the mass spectrum in the quantum YM theory with periodic boundary con- 
ditions ]^0| . Another possible direction is a clarification of a connection of the present 
approach with the quasiclassical quantization of periodic orbits which hopefully could be 
useful for the quantum chaos studies. In any case, a lot of work is needed before the 
suggested method proves or disproves its power. 
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Note added When this paper has been completed, I became aware of the recent paper 
[31] where similar ideas were developed and applied to the quantum anharmonic oscillator. 
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